Measurement induced chaos with entangled states 
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The dynamics of an ensemble of identically prepared two-qubit systems is investigated which is 
subjected to the iteratively applied measurements and conditional selection of a typical entanglement 
purification protocol. It is shown that the resulting measurement-induced non-linear dynamics of 
the two-qubit state exhibits strong sensitivity to initial conditions and also true chaos. For a special 
class of initially prepared two-qubit states two types of islands characterize the asymptotic limit. 
They correspond to a separable and a maximally entangled two-qubit state, respectively, and their 
boundaries form fractal-like structures. In the presence of incoherent noise an additional stable 
asymptotic cycle appears. 
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Introduction — Entanglement is at the heart of quan- 
tum physics since its discovery [l|. However, it was only 
recently that the focus has been put on entanglement 
as a resource for quantum communication and quantum 
information processing Various protocols have been 
developed to detect Q , generate and distill entanglement 
From an ensemble of identical quantum states, one 
can produce an ensemble yielding higher degree of entan- 
glement by unitary transformations, measurements and 
selection conditioned on measurement outcomes. These 
non-linear processes are referred to as entanglement dis- 
tillation or purification. Entanglement distillation proto- 
cols play a crucial role in increasing the quality of com- 
munication channels and have also been used to define 
the degree of entanglement in an operational sense. 

The phenomenon that sensitivity to initial conditions 
leads to chaotic dynamics in classical physics is well- 
known. Similar phenomena in closed quantum systems 
are, however, excluded by the quantum unitary evolution 
Q . In the case of open quantum systems the restriction 
to unitarity is lifted. The evolution of an open quantum 
system is sometimes pictured as additionally having an 
environment that performs generalized measurements on 
it. In general, any type of measurement makes the evo- 
lution non-unitary. Entanglement purification, i.e. se- 
lection of certain systems from an ensemble of identical 
systems based on the results of partial measurements, can 
be regarded as a generalized feedback mechanism. The 
non-linear dynamics resulting from such generalized feed- 
back has been shown to lead in certain cases to a chaos 
0. This type of chaos is essentially different from that 
arising from the stochastic dynamics of a continuously 
measured open quantum system Q, which can become 
chaotic in the semi-classical regirne while still showing 
signatures of quantum behaviour [8| . 



The generalized feedback resulting from measurement 
based selection plays a crucial role in the case of vari- 
ous entanglement purification protocols 41 . Applying the 
purification protocol 9] on an ensemble of single qubits 
prepared in identical pure states, after each iteration step 
the remaining, selected ensemble of qubits will again be 
in identical, pure quantum states. The dynamical evo- 
lution can be characterized by a rational non-linear map 
over the extended complex plane representing the pure 
states. This map has been proven to lead to truly chaotic 
behaviour [6] . A single qubit is the simplest quantum sys- 
tem, and also lacks all genuine non-local quantum proper- 
ties. Thus the above mentioned study left a fundamental 
question open, which in turn raises naturally in the con- 
text of entanglement purification: Can we find sensitiv- 
ity to initial conditions for genuine mulipartite quantum 
properties, in particular for entanglement, when only lo- 
cal operations and classical communication (LOCC) are 
applied? 

In this Letter we focus on a particular entanglement 
purification protocol Q, and demonstrate the existence 
of true chaos which manifests also in the evolution of 
entanglement. An important feature of the protocol is 
that it maps pure states onto pure states, moreover it 
may also increase purity of initially mixed states. Our 
analysis is accomplished by showing that the convergence 
to a fully entangled or a separable asymptotic attractor 
can be sensitive to the initial state. 

Measurement based non-linear dynamics — The pure 
quantum state state of a system consisting of a pair of 
qubits can be expressed in the computational basis as 
I*) = ci|00) + C2IOI) + C3IIO) + C4III). We consider 
an ensemble of qubit pairs prepared in the state jvP) as 
the input to the entanglement purification protocol Q 
schematically depicted on Fig. [1] This protocol realizes 
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FIG. 1. Schematic representation of one iteration step of the 
entanglement purification process. The diamond shaped ele- 
ments denote acceptance of the pair of qubit only when both 
measurements yield 0. The depicted procedure is used to pre- 
pare from an ensemble E made up of pairs of qubits all in a 
state IV'), an ensemble E' made up of pairs of qubits in the 
state I t/)' ) . During the procedure half of the pairs in E are com- 
pletely used up, while a portion of the other half is retained 
depending on the measurement outcome. Here we depicted 
the special case U — H (g) H , where H is the Hadamard gate. 



a non-linear transformation of the quantum state | VP) ac- 
cording to Ci Ncf , where is a necessary normaliza- 
tion factor. The non-linearity is due to the generalized 
feedback realized by the measurement-based conditional 
selection. By following each non-linear transformation 
by a unitary U, further non-trivial dynamics can be gen- 
erated. A complete iteration can be expressed as the 
transformation 



(1) 



Note, that only local unitary transformations comply 
with the concept of entanglement purification. For the 
following analysis we fix the unitary transformation to 
he U — H ® H , where iJ is a one qubit Hadamard gate 

For Hilbert spaces of dimensionalities of more than 
two, such as for that of a two qubit system, the represen- 
tation of a pure state requires a vector of several com- 
plex numbers, and the non-linear dynamics is described 
by a non-linear map on this higher dimensional complex 
space. The mathematics of non-linear dynamical maps in 
several complex variables is substantially more involved 
and much less understood than the same for a single com- 
plex variable. Even the existence of chaotic regions is a 
nontrivial question fl^. 

We are interested in the evolution of genuine multi- 
partite properties of the system under the iterative non- 
linear dynamics. Since this protocol describes an entan- 
glement purification protocol, of all such properties, en- 
tanglement is of the most concern. Thus we consider 
states parametrized by the complex number including in- 
finity C (C g C) of the form 



where the normalisation factor reads M{C) = (1 + 
|(^|2^-i/2_ The degree of entanglement of this state is 
completely determined by C via the binary entropy ivaic- 



m. 



By studying the evolution of these states we can learn 
about the general multi-partite properties of this proto- 
col. The analytic treatment of this evolution is greatly 
simplified by the fact that the states of the form of Eq. ([2]) 
are invariant under two successive iterations, in particu- 
lar yielding 
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Thus the description of the dynamics for this class of 
initial states simplifies to a non-linear map of a single 
complex variable. The function g{C,) : C i—)- C is a fourth 
order rational function generating a one variable com- 
plex dynamical map on the Riemann sphere. The map 
is chaotic in the sense that the corresponding Julia set is 
non- vacuous [l^ . Fourth order maps in general can lead 
to rather involved behaviour. We base our analysis on 
the observation that the fourth order rational map can 
be written as a composition of a second order rational 
function with itself 5(C) = /(/(C)) = /°^(C), where 



/(C) = 



(5) 



Moreover, the same function also describes the quantum 
state after an odd number of iterations, according to 



^ M (|00) + 111) + /°'"+\C)(|oi) 



|io))). 

(6) 

Thus the iterative dynamics of initial states of the form 
Eq. (12) is equivalent to the iterative dynamics generated 
by the second order rational map defined in Eq. ([5|). In 
the following we will examine this mapping in some de- 
tail. 

Stable cycles of the dynamics — The long_term be- 
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haviour of a rational map can be analyzed 12| by fol- 
lowing the orbits of its critical points (points where the 
derivative of the map vanishes /'(Cc) — 0). The second 
order map in Eq. ([5|) has two critical ponts: Cci = and 
Cc2 = 00. The first one, Cci is part of the superattracting 
cycle {0, 1}, while the second one lands on the same cycle 
after two iterations. Therefore, the map has one stable 
(superattractive) cycle. Translated back to the language 
of states, we have to take into account that the function 
/ describes a state in the same form only in every sec- 
ond step. Thus, depending on the parity of steps when 
reaching the first element of the cycle, we have two dis- 
tinct cases. 



3 



In the first case, when /°^"(C) = 1, the corresponding 
state reads 

|vI/(2«)) = (100) + 111)) = !<!>+). (7) 

After an even number of steps, one reaches a fully entan- 
gled state, the Bell state 1$+). Since this state is invari- 
ant under both the non-linear transformation S and the 
unitary operation H<g)H,we have |vI/(2")) = |\i/(2n-(-i)^ ^ 
|$+). Therefore, the Bell state \^~^) is an asymptotically 
stable fixed point of the dynamics. 

In the second case, when /°^"(C) — 0, the correspond- 
ing state reads 



1^(2")) = |00) 



(8) 



The above product state is not invariant under the iter- 
ative dynamics, however, any subsequent step will leave 
the state completely separable. In particular, after an 
odd number of steps we find 



I* 



■2n+l) 



|1))®(|0) + |1)) 



(9) 



Therefore, the second stable cycle is of length two, both 
of its members are separable pure states. 

Due to a theorem on rational maps 1^, the degree of 
the rational function determines the maximum number 
of stable cycles. For a rational map of degree two, at 
most two stable cycles can exist. In our case, the single 
stable cycle of the rational map / can lead to two different 
stable cycles of the dynamics, depending on the parity of 
the number of steps when approaching the first element 
of the cycle. Thus, we have found all possible stable 
cycles of the dynamics restricted to initial pure states of 
the form 

Sensitivity to initial states — The two possible stable 
cycles of the dynamics are very different. One of them 
is a single, completely entangled pure state, a Bell state, 
while the other is an oscillation between two separable 
pure states of the two qubits. We ask now the question, 
what are the initial states converging to each of the stable 
cycles. 

Let us first discuss the case of real values for the pa- 
rameter The function / maps real numbers to real 
numbers, thus it can be restricted to M. The members of 
the stable cycle {0, 1} are also real numbers. We can now 
determine the basin of attraction for the two cases of con- 
vergence, i.e. convergence to after an even or an odd 
number of steps, which we shall call even-zero or odd- 
zero convergence, respectively. The immediate neigh- 
bourhood of the fixed point 0, belonging to even-zero 
convergence, is determined by the equation I /°^((^) I < |(|, 
with the condition \Q < 1. The corresponding equation 
can be explicitly solved yielding |^| < (a where 



The preimages of the interval {—(a, Ca) belong to an odd- 
zero convergence region. By solving the corresponding 
equation we find two distinct intervals of odd-zero con- 
vergence Ca < ICI < C-B where 
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It is easy to see that the region |C| > Cs is mapped after 
one iteration to the region (a < |CI < Cb, thus it belongs 
to even-zero convergence. To summarize the behaviour 
of the map restricted to the reals, we have found that re- 
gions of odd-zero and even-zero convergence follow each 
other, these open sets belong to the Fatou set. The bor- 
der point (^A is a repulsive fixed point of the map, while 
the other border points {—Cat—CbXb} are preimages of 
(a- These four points belong to the Julia set. 

A similar analysis can be repeated for the map / with 
domain C. The preimages of with their open small 
neighbourhoods provide regions of odd-zero or even-zero 
convergence, forming the Fatou set. The Fatou set is an 
open set, while the complementary closed set is the Ju- 
lia set. A small neighbourhood of the origin will belong 
to even-zero convergence, a sufficient condition for this 
is that |/(/(C))| < ICI- From this condition we can de- 
termine the maximum radius Cc of a circle around the 
origin belonging to even-zero convergence. The radius Cc 
can be calculated explicitly from an algebraic equation, 
its numerical value is Cc ~ 0.475. The first order pre- 
images of are 1 and —1. Thus, 1 and —1 together with 
a region around them belong to odd-zero convergence. 
The preimage of the circle of convergence around zero 
determines an immediate region of convergence around 1 
and —1. In a similar manner one can continue this pro- 
cess and determine the next order preimages and regions 
of convergence around them. The numerically calculated 
convergence regions are shown in Fig. [2] We can clearly 
recognize a chessboard-like structure of green and blue 
islands, belonging to convergence to fully entangled or 
completely separable states, respectively. The islands fol- 
low a self similar structure with decreasing size. A Julia 
set is formed by the border points between the two colors, 
with no stable cycles in it. The asymptotic entanglement 
of the two-qubit system behaves chaotically, it is sensi- 
tive to the initial state on arbitrary small scales. Since 
this is an asymptotic map, it involves an infinite-in-time 
limit process. 

Mixed initial states — Up to now, we have considered 
pure initial states. Adding noise to the initial state, de- 
scribed by a density operator, will alter the dynamics. 
We will test the sensitivity to a certain type of noise by 
adding the unit matrix to the density matrix representing 
the initial pure state 



p(C,A)-A|*(C))(*(C)| 



1- A 



-1, 



(12) 



Ca = (a ~ 1 ~ 2/a)/3 with a = (17 + 3\/33)^^^ . (10) where A G [0, 1], and 1 stands for the unit operator acting 
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FIG. 2. Convergence of initial states of Eq. ((2|) to the two 
limiting cycles as a function of the complex parameter (. The 
blue color denotes convergence to the maximally entangled 
state of Eq. (O, and the green color the convergence to the 
cycle of the two separable states in Eqs. ^ and @. 



on the Hilbert space of the two qubits. Since the dynam- 
ics is no longer restricted to pure states, we can expect 
that further stable fixed cycles will appear, containing 
mixed states. The stability of the fixed cycles can be 
proven in any convenient representation. We chose the 
Fano representation where the real expansion coefficients 
with respect to the 16 generalized Pauli matrices repre- 
sent an arbitrary density matrix 13]. Moreover, the Fano 
representation is convenient for numerical simulation of 
the dynamics also. By calculating the eigenvalues of the 
Jacobian in the Fano representation for each numerically 
found cycle, we concluded that among them the only sta- 
ble cycle is the length two cycle {pi,p2}, where 

pi = i(|00)(00| + |ll>(ll|), (13) 

P2^\ {mm + |ii>(ii| + (101) + |io)) ((Oil + (lOD) . 

The same calculation for the cycles known from the pure 
initial state case indicated their stability against pertur- 
bation by arbitrary mixed states. 

In Fig. |3]we show the convergence towards the stable 
cycles of the mixed state dynamics. The third, mixed 
stable cycle {pi,p2} denoted by yellow color washes out 
the fine structure of the pure state picture. Islands of 
purification towards both fully entangled and completely 
separable states remain visible even for the value A = 0.75 
of the mixing parameter. Understanding the structure of 
convergent areas requires further studies. An interesting 
question is, whether the number of the islands purify- 
ing to entangled states is finite or it possesses a fractal 
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FIG. 3. Convergence of mixed initial states p(A, of Eq. (|12p 
to the three limiting cycles as a function of the complex pa- 
rameter (, at A = 0.75. The blue and green colours denote 
convergence to the same cycles as on Fig. [21 while the yel- 
low colour denotes convergence to the separable mixed cycle 
in Eq. (|13l) . As A is increased from to 1, the correspond- 
ing convergence plots gradually turn from that on Fig. [2] to 
uniform yellow. 



structure. 

Conclusions — We have demonstrated that entangle- 
ment in a quantum system can evolve truly chaotically, 
exhibiting sensitivity to the initial condition. Our re- 
sults give an insight into the properties of the general 
pure state dynamics of a protocol that is described by 
a non-linear dynamical map in three complex variables. 
We have found that extending the space of initial states 
to a special class of mixed states, a new, mixed attracting 
cycle appears. It would require further studies to decide, 
whether adding other types of noise to the initial state 
would introduce additional attractors. 
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